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Abstract. The formation of H2 and HD molecules on interstellar dust grains is studied using 
rate equation and master equation models. Rate equations are used in the analysis of laboratory 
experiments which examine the formation of molecular hydrogen on astrophysically relevant 
surfaces. However, under interstellar conditions, rate equations are not suitable for the calcu- 
lation of reaction rates on dust-grain surfaces. Due to the low flux and the sub-micron size of 
the grains, the populations of H and D atoms on a single grain are likely to be small. In this 
case the reaction rates are dominated by fluctuations and should be calculated using stochastic 
CIh 1 methods. The rate of molecular hydrogen formation in interstellar clouds is evaluated using the 
D ' master equation, taking into account the distribution of grain sizes. 
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1. Introduction 

The formation of molecular hydrogen in the interstellar medium is a process of fun- 
t ■ damental importance in astrophysics. H2 molecules contribute to the cooling of clouds 
\ during gravitational collapse and participate in reactions that produce more complex 
■ molecules ( |Duley fc Williams 1984||Wilriam~ 998). H2 cannot form in the gas phase ef- 
l/-) \ ficiently enough to account for its observed abundance ( |Gould fc Salpeter 1963| ). It was 
<^ ■ thus proposed that H2 formation takes place on dust grains that act as catalysts. The 
**^h \ process of H2 formation on grains can be broken up into several steps as follows. An H 
O atom approaching the surface of a grain has a probability £ to stick and become adsorbed. 
' I The adsorbed H atom (adatom) resides on the surface for an average time in (residence 
;| ■ time) before it desorbs. In the Langmuir-Hinshelwood mechanism, the adsorbed H atoms 
I quickly equilibrate and diffuse on the surface of the grain either by thermal activation or 
d • by tunneling. When two adsorbed H atoms encounter each other, an H2 molecule may 

form HWilliams 1968|Hollenbach fc Salpeter 1970|Hollenbach fc Salpeter 1971IHollenbach et al. 1971ISmoluchowski 1981|Arono-w 
The steady state production rate of molecular hydrogen, Rn 2 (cm -3 s _1 ) can be expressed 

by 



i?H 2 = ^naVHajng, (1.1) 

where 71h (cm -3 ) and «h (cm s _1 ) are the number density and the speed of H atoms 
in the gas phase, respectively, a (cm 2 ) is the average cross-sectional area of a grain and 
n g (cm -3 ) is the number density of dust grains. The parameter 7 is the fraction of H 
atoms striking the grain that eventually form a molecule, namely 7 — £77, where 77 is 
the probability that an H adatom on the surface will recombine with another H atom 
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to form H2. The probability £ for an H atom to become adsorbed on a grain surface 
covered by an ice mantle was calculated by Buch and Zhang (1991) and Masuda et al. 
(1998). They found that £ depends on the surface temperature and on the energy of the 
irradiation beam. For a surface at 10 K and beam temperature of 350 K, Masuda et al. 

(1998) obtained a sticking coefficient around 0.5. 

A simplified version of Eq. I|l.lfl is often used for the evaluation of the H2 formation 
rate in models of interstellar chemistry. It takes the form 

i?H 2 = R ■ nun, (1.2) 

where R (cm 3 s _1 ) is the rate coefficient and n = ?ih + 2nn 2 (cm' 3 ) is the total density 
of hydrogen atoms in both atomic and molecular (nn 2 ) form. In previous studies it was 
assumed that the total grain mass is about 1% of the hydrogen mass in the cloud, with 
a single grain size of r = 0.17/zm and a grain-mass density of p g — 2 (gram cm~ 3 ) 
( Hollc nbach et al. 1971J1 . Under these assumptions, for gas temperature of 100 K and 
7 = 0.3, one obtains R ~ 10~ 17 (cm 3 s _1 ), in agreement with observations for diffuse 
clouds IjJura 1975(1 . More recent observations in photon dominated regions, indicate that 
in these regions the rate coefficient should be R ~ 10~ 16 (cm 3 s _1 ), in order to account 
for the observed H 2 abundance CH abart et al. 2003jl . 

2. Laboratory Experiments 

In the last few years, the formation of molecular hydrogen on interstellar dust analogues 
was studied in a series of laboratory experiments IjPirronello et al. 1997alPirronello et al. 1997blPirronello et al. 1999IRoser et al. 
In these experiments, beams of H and D atoms were irradiated on the surface. The 
production of HD molecules was measured both during irradiation and during a sub- 
sequent temperature programmed desorption (TPD) experiment. In this experiment, 
the temperature of the sample is raised quickly to either desorb particles (atoms and 
molecules) that got trapped on the surface or to enhance their diffusion and subsequent 
reaction and desorption. The formation of molecular hydrogen was studied on samples of 
polycrystalline olivine HPirronello et al. 1997alPirronello et al. 1997b|) . amorphous car- 
bon IjPirronello et al. T 999 ) and various types of amorphous ice samples, namely high 
density ice (HDI) and low density ice (LDI) l|Manico et al. 2001IRoser et al. 2002IHornekaer et al.2003"|l . 
The results of the TPD experiments, particularly on polycrystalline olivine and HDI, 
show second order kinetics for low coverage of H and D atoms on the surface. In second 
order kinetics the peak shifts towards higher temperature as the coverage is reduced. 
This is due to the larger distancees that atoms need to diffuse before they encounter 
each other and recombine. This indicates that on these samples the surface mobility of 
H and D atoms is dominated by thermal hopping rather than tunneling. If tunneling 
was dominant, molecules could form efficiently during irradiation, resulting in first order 
kinetics. 

The results of the TPD experiments were analyzed using rate equation models IjKatz et al. 1999ICazaux & Tielens 2002ICazaux 
The analysis involves fitting the parameters of the rate equations to the experimental 
TPD curves. The most important parameters are the activation energy barriers for dif- 
fusion and desorption of hydrogen atoms on the surface. The models used by Katz et al. 

(1999) and Perets et al. (2005) take into account only physisorption sites and thermal 
hopping between them. The models of Cazaux and Tielens (2002;2004) also take into 
account tunneling and chemisorption sites. However, the two analyses are in agreement 
on the fact that tunneling between physisorption sites appears to be too slow to play an 
important role in the TPD experiments (Vid ali et al. 2005}) . Furthermore, recent experi- 



The formation of H2 and HD 



3 



ments and calculations indicate that there are energy barriers which are likely to prevent 
adsorbed H and D atoms from entering chemisorption sites at these low temperatures. 

Using the parameters obtained in the laboratory experiments, the rate equations en- 
abled the extrapolation from laboratory time scales to astrophysical time scales. Using 
rate equations, the efficiency of molecular hydrogen formation on these surfaces under 
interstellar conditions and its dependence on the flux and grain temperature were cal- 
culated. It was found that the recombination efficiency strongly depends on the grain 
temperature. Each sample exhibits a narrow window of high efficiency along the temper- 
ature axis. For the polycrystalline olivine sample, this window is between 7 — 9K, while 
for amorphous carbon it is between 11 — 17 K (Vida li et al. 20 051. 



3. The Rate-Equation Model 

Consider a diffuse interstellar cloud with a density njj (cm -3 ) of H atoms, and a density 
rig of dust grains. The typical velocity vjj (cm s _1 ) of H atoms in the gas phase is given 

by 



/ 8 k B T gas 

V H = \ — ! I 3 - 1 ) 

V 7T TOR 

where mjj = 1.67- 10~ 24 (gram) is the mass of an H atom and T gas is the gas temperature. 
To evaluate the flux of atoms onto grain surfaces we will assume, for simplicity, that the 
grains are spherical with a radius r (cm). The cross section of a grain is a = irr 2 and its 
surface area is 4cr. The density of adsorption sites on the surface is s (sites cm -2 ) and 
the number of adsorption sites on a grain is S = 4irr 2 s. The flux F H (atoms s _1 ) of H 
atoms onto the surface of a single grain is given by Fh = n^v-aa. 

H atoms that collide and stick to the surface hop as random walkers between adjacent 
sites until they either desorb or recombine into molecules. The desorption rate of an H 
atom is 



W H = u-exp(-E 1 /k B T), (3.2) 

where v is the attempt rate (typically assumed to be 10 12 s _1 ), E\ is the activation energy 
barrier for desorption of an H atom and T is the surface temperature. The hopping rate 
of an H atoms is 



an = v ■ exp(-E /k B T), (3.3) 

where Eq is the activation energy barrier for H diffusion. Throughout this paper we use 
the parameters obtained experimentally for amorphous carbon, namely the activation 
energies are Eq = 44.0 meV and E\ — 56.7 meV fKatz et al. 19 99). and the density of 
adsorption sites on the surface is s ~ 5 x 10 13 (sites cm -2 ) IjBiham et al. 200"T)l . For the 
density of hydrogen atoms in the gas phase we take tt-h = 10 (atoms cm -3 ). The gas 
temperature is taken as T gas = 90 K, thus v B = 1-37 x 10 5 (cm s -1 ). These are typical 
values for diffuse interstellar clouds. 

The number of H atoms on the grain is denoted by TVh and its expectation value 
is denoted by (Njj). The rate An — an/S is approximately the inverse of the time t s 
required for an atom to visit nearly all the adsorption sites on the grain surface. This is 
due to the fact that in two dimensions the number of distinct sites visited by a random 
walker is linearly proportional to the number of steps, up to a logarithmic correction 
(jMontroll fc Weiss 1965"jl . 
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The formation of H2 molecules on interstellar dust grains can be described by the rate 
equation 

^^=F ii -W ii {N H )-2A H (N ii ) 2 . (3.4) 

The first term on the right hand side describes the flux of H atoms, the second term 
describes the desorption of H atoms and the third term describes the diffusion and 
recombination. Here we assume, for simplicity, that all H2 molecules desorb from the 
surface upon formation. The production rate i?n 2 (cm -3 s -1 ) of H2 molecules is given 

by 

i?H 2 =A H {N u } 2 n g . (3.5) 
Solving the rate equations under steady state conditions for astrophysically relevant 
values of the flux and temperature, it was found (Biham & Lipshtat 2002) that molecular 
hydrogen formation is efficient only within a narrow window of grain temperatures, To < 
T < T\, where 

T ° = k^sJK) (3 ' 6) 

and 

Tl = k B \n{vS/F^y (3J) 
At temperatures below T the mobility of H atoms on the surface is very low, sharply 
reducing the production rate. At temperatures above T% most atoms quickly desorb before 
they encounter each other and form molecules. 

Rate equations are an ideal tool for the simulation of surface reactions, due to their 
simplicity and high computational efficiency. In particular, they account correctly for the 
temperature dependence of the reaction rates. However, in the limit of small grains under 
low flux they become unsuitable. This is because they ignore the fluctuations as well as the 

discrete nature of the populations of atoms on the grain ( |Charnley et al. 1997ICaselli et al~ 998 Shalabie a et al. 1998IStantcheva ■ 
For example, as the number of H atoms on a grain fluctuates in the range of 0, 1 or 2, 
the H2 formation rate cannot be obtained from the average number alone. This can be 
easily understood, since the recombination process requires at least two H atoms simul- 
taneously on the surface. In this limit, the master equation approach is needed in order 
to evaluate the H 2 formation rate. 



4. The Master Equation 

Consider a grain exposed to a flux of H atoms. The probability that there are Nh 
hydrogen atoms on its surface is given by Ph(-^Vh), where = 0, 1, 2, . . . , S. The master 
equation provides the time derivatives of these probabilities. It is a set of coupled linear 
differential equations of the form 

Ph(N u ) = F H [P H (N H - 1) - Jfc(JV H )] + VK H [(JV h + l)PH(iVH + 1) - N H P K (N H )} 

+ A H [(JV h + 2)(JV H + 1)P K (N H + 2) - jV H (JV H - 1)Ph(/V H )] . (4.1) 

Each equation includes three terms. The first term describes the effect of the incoming 
flux _Fh. The probability Ph(JVh) increases when an H atom is adsorbed on a grain that 
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already has Ah — 1 adsorbed atoms [at a rate of Fh -Ph (Ah — 1)] , and decreases when a new 
atom is adsorbed on a grain with Ah atoms on it [at a rate of FhPh(Ah)]. The second 
term includes the effect of desorption. An H atom that is desorbed from a grain with 
A H atoms, decreases the probability Ph(A h ) [at a rate of AhIFhPr^Ah)], and increases 
the probability Ph(Ah — 1) at the same rate. The third term describes the effect of H2 
formation. The production of one molecule reduces the number of H atoms on the surface 
from Ah to Ah — 2. For one pair of H atoms the recombination rate is proportional to 
the sweeping rate multiplied by 2 since both atoms are mobile simultaneously. This 
rate is multiplied by the number of possible pairs of atoms, namely Ah (Ah — l)/2. Note 
that the equations for Ph(0) and Ph(1) do not include all the terms, because at least 
one H atom is required for desorption to occur and at least two for recombination. The 
rate of formation of H2 molecules on a single grain, -Rfj™ 11 (cm -3 s -1 ), is thus given by 

R^ = A n [(N*)-(N H )] (4.2) 

where 



<JVh> = £ ^h(Ah). (4.3) 

N H =0 

is the fcth moment of the distribution Ph(Ah). 

The time dependence of i? H ' 2 am can be obtained by numerically integrating Eqs. I|4.1(l us- 
ing a standard Runge-Kutta stepper. For the case of steady state, namely Ph(A h ) = for 
all A H , an analytical solution for Ph(A h ) was obtained, in terms of Ah/Wh and W-h/Fh 
IjGreen et al. 2 001 Biham & Lipsht at 2002| ). The rate of formation of H2 molecules on a 
single grain with S adsorption sites is given by QGreen et al. 2001|Biham fc Lipshtat 2 002 ) 



i?!™ 11 = -ii j ' , (4.4) 



where I x {y) is the modified Bessel function. 

In Fig. 1 we present the H2 production rate on a single grain of amorphous carbon 
vs. the grain size, which is parametrized by the number of adsorption sites, S. The 
conditions, specified in Sec. 3, are typical in diffuse clouds and the grain temperature 
is T = 18 K. It is observed that for large grains the rate equation results coincide with 
those of the master equation. However, for small grains the rate equation over-estimates 
the production. The master equation accounts correctly for the rate of H2 formation on 
grains of all sizes. 



5. Effects of the Grain Size Distribution 

Interstellar dust grains exhibit a broad size distribution between several nano-meters 
and fractions of a micron. This distribution can be approximated by a power law of the 
form IjMathis et al. 1977|Weingartner fc Draine 2001| ) 

n g (r) = cr- a (5.1) 

where n g (r)dr (cm -3 ) is the density of grains with radii in the range (r, r + dr) in the 
cloud. This distribution is bounded between the upper cutoff r max = 0.25/im and the 
lower cutoff r m i n = 5nm. The exponent is a ~ 3.5 IjDraine fc Lee 1984(1 . Clearly, this 
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Figure 1. The production rate of H2 molecules on a single grain of amorphous carbon with 
S adsorption sites. For large grains the rate equation results (dashed line) are in agreement 
with the master equation (solid line) , while for small grains the rate equation over-estimates the 
production. The parameters are specified in Sec. 3. 



distribution is dominated by the small grains. This fact has important implications on 
the formation rate of molecular hydrogen. The total surface area is much larger than 
in the case of a homogeneous population of grains with radii of the order of 0.25/im, 
with the same total mass density of about 1% of the hydrogen mass density. This tends 
to enhance the H2 formation rate. However, below some grain size the efficiency of H2 
formation is sharply reduced. Due to these competing effects the production of molecular 
hydrogen should be integrated over the whole range of grain sizes. 

In Fig. 2 we present the rate coefficient R (cm 3 s -1 ) vs. grain temperature for the 
case of a power-law distribution of grain sizes (solid line) and for the case in which all 
grains are of radius r = 0.17/xm, with the same total mass of grains. It turns out that 
within the temperature window in which H2 formation is efficient, the increase in the 
total surface area due to the broad size distribution, enhances the production rate by 
about an order of magnitude compared to a homogeneous population of sub-micron size 
grains flLipshtat fc Biham 2QQ5} . This is in spite of the reduced efficiency of H2 formation 
on small grains. 
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Figure 2. The rate coefficient R (cm 3 s _1 ) of Eh formation vs. grain temperature, obtained 
from the master equation for the case of a power law distribution of grain sizes (solid line) . The 
dashed line shows the rate coefficient for the case in which the same mass of grains consists of 
a homogeneous population of "classical grains" of radius r = 0.17/im. 

6. Effects of Grain Porosity 

The laboratory experiments on interstellar dust analogues have shown that H 2 for- 
mation on dust grains is efficient within a range of grain temperatures below 20 K. 
However, high abundances of H2 have also been observed in warmer clouds, including 
photon-dominated regions (PDRs), where grain temperatures may reach 50 K. The pos- 
sible effect of grain porosity on H2 formation was recently examined using a rate equation 
model (|Perets fc Bi ham 2005). It is found that porosity extends the temperature range 
in which H2 formation is efficient towards higher temperatures. This is because H atoms 
that desorb from the internal surfaces of the pores may re-adsorb many times and thus 
stay longer on the surface. However, this extension is expected to be of only several 
degrees (K) and by itself cannot account for the formation of molecular hydrogen in 
PDRs. 

Another effect of porosity that may be relevant to molecular hydrogen formation is that 
porous grains may be significantly colder than predicted by models which use compact, 
spherical grains l|Blanco fc Bussoletti 19 80 Voshchinniko v"et al. 19 99) . While models of 
spherical grains predict grain temperatures in the range of 30 — 50 K in PDRs, the actual 
temperatures of porous interstellar grains may thus be significantly lower. In the most 
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favorable cases, the grain temperatures may go down to the range of 15 — 30 K for large 
grains and to 10 — 20 K for small grains IjBlanco fc Bussoletti 1980p . 

7. HD Formation on Grains 

Unlike the H2 molecules in the interstellar medium which do not form efficiently in the 
gas phase, HD molecules have an efficient gas-phase formation mechanism. It consists of 
a sequence of reactions in which an ionized D atom reacts with an H2 molecule to form 
HD. The possibility that in diffuse clouds some HD molecules are also formed on grain 
surfaces was recently examined |Lipshtat et al. 2004| ). Due to the lower abundance of D 
atoms vs. H atoms (by a factor of about 10 _b ), the master equation is required for such 
calculations even on relatively large grains. In case that D atoms bound more strongly to 
the surface than H atoms, the temperature window in which HD formation is efficient is 
extended towards higher temperatures, compared to the range of efficient H2 formation. 
The assumption that D atoms are more strongly bound to the grains is plausible, but so 
far has not been confirmed experimentally. 

8. Complex Reaction Networks of Multiple Species 

The master equation can be generalized to describe complex reaction networks of multi- 
ple species on grain surfaces. Recently, it was applied to the study of the reaction network 
that lead to the formation of H2CO and CH3OH on grains (Stantcheva et al. 2002 Stantchcva & Herbst 2003). 
Such surface reaction networks are dominated by H atoms because hydrogen is the most 
abundant specie and also because it reacts with most of the other atomic and molecu- 
lar species adsorbed on grain surfaces. Furthermore, hydrogen atoms hop on the surface 
much faster than other adsorbed species and thus quickly react with other species even 
at very low coverages. For such reaction networks, the dynamical variables in the master 
equation are the probabilities P(Nh, No, Noh, ■ ■ ■) that there are Nx copies of specie 
X on the grain. Setting suitable cutoffs, each population of each specie is limited to 
iVx = 0, 1, . . . , Af™ ax . The problem is that as the number of reactive species increases, the 
number of coupled equations in the master equation increases exponentially, limiting the 
feasibility of this approach. Using the fact that the reaction networks are usually sparse, 
it is possible to dramatically reduce the number of equations ( Lipshtat et al. 2004 ) . This 
makes the master equation applicable even for complex surface reaction networks with a 
large number of reactive species. 

The master equation can be solved either by direct integration, as described above, 
or by Monte Carlo (MC) methods ( |Charnley 2001| ). A significant advantage of direct 
integration is that the equations can be easily coupled to the rate equations of gas-phase 
chemistry. Furthermore, a single run of the direct integration method provides the entire 
probability distribution of populations of adsorbed species for the simulated conditions 
as well as the reaction rates. This is in contrast with the MC methods, which typically 
require large computational resources due to the need to accumulate much statistical 
information. 



9. Summary 

We have presented the master equation approach for the evaluation of molecular hy- 
drogen formation rates on interstellar dust grains. This approach is required in the limit 
of small grains and low flux where the population of hydrogen atoms on each grain is 
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small. In this limit, the reaction rate is dominated by fluctuations. Thus, the rate equa- 
tion approach, which incorporates the mean-field approximation fails. Applications of the 
master equation approach for the evaluation of the formation rates of H 2 , HD as well as 
more complex molecules are discussed. 
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Discussion 

Hornekaer: Will the inclusion of chemisorption sites on dust grains turn even small 
interstellar grains into efficient contributors to H2 formation in photon dominated regions 
(PDRs) due to the fact that thermal spikes could activate recombination of chemisorbed 
H atoms? 

Biham: In our analysis so far we have not included fluctuations in grain temperatures. 
Such fluctuations are expected to be significant only for grains of radii smaller than about 
10 ran. Chemisorbed H atoms are expected to contribute to H2 formation when grain 
temperatures reach several hundred degrees (K). However, a careful analysis is required 
in order to quantify the contribution of such mechanism. 

AlKAWA: You concluded that H migrates by thrmal hopping, while Cazaux fits the same 
experimental data assuming tunneling. Would you tell me why you need to exclude 
tunneling and why you and Cazaux reached different conclusions on H-migration? 

Biham: The temperature programmed desorption (TPD) curves obtained for the poly- 
crystalline olivine sample (Katz et al. (1999) exhibit second order kinetics at lower cov- 
erage, namely, as the coverage decreases, the TPD peak shifts to higher temperatures. 
This indicates that the HD molecules are formed by thermal hopping during the TPD 
run. If tunneling alone provided sufficient mobility for H atoms to diffuse and recombine 
during irradiation, when the surface temperature is around 6 K, all TPD curves would 
be of first order kinetics. For other samples, such as the amorphous carbon, the evidence 
for second order kinetics is not as strong and tunneling is not ruled out completely. 



